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Abstract 



The first in m/M and fourth in Za (pure recoil) correction to a hydrogen 
energy level is found. This correction comprises two contributions, one coming 
from the atomic scale, the other from distances of the Compton wavelength order. 
Two different perturbation schemes are used to calculate the former. One of them 



exploits as unperturbed the solution to the Dirac-Coulomb problem, the nucleus' 
slow motion being the source of the perturbation. The alternative scheme treats 
both the electron and the nucleus as slow, while relativistic effects are consid- 
■ ered perturbatively. The short-distance contribution is found in the Feynman and 

Coulomb gauges. Recent results for P levels are confirmed, in contrast with those 
for S levels. Numerically, the shift equals 2.77 kHz for the ground state and 0.51 
kHz for 2S state. 



1 Introduction 



To determine the proton charge radius with a percent accuracy from the value of the 
hydrogen Lamb shift, the latter should be known both experimentally and theoretically 
with a precision of one kHz. Recently completed calculation of the order ma 2 (Za) 5 
corrections ]I[ leaves among the effects of possible phenomenological interest the pure 
recoil correction of the order m 2 (Za) 6 /M arising due to interference between the nucleus' 
recoil and relativistic effects in the motion of the electron. The present paper is devoted 
to the calculation of this correction for an arbitrary state of the hydrogen atom. 

Recently this correction for P states was found ||. In those states, as well as in all 
states with nonzero angular momenta, the correction proves to be saturated by a con- 
tribution coming from the atomic scale. Hence one can use there the standard quantum 
mechanical perturbation theory for the effective operators describing relativistic effects. 
Matrix elements of the effective operators arising in the perturbation theory converge 
at small distances thus testifying a posteriori that the used 'nonrelativistic' approach is 
correct at the given order of a for states with nonzero I. 

An attempt to apply the same approach to S states, whose wavefunctions do not 
vanish at the origin, leads to matrix elements diverging at small distances. In fact, among 
the effective operators one finds those depending on r as r -3 and even r -4 ||. As for 
the latter, for S states the operator r~ 4 is equivalent (modulo a nonsingular operator) to 
the sum of operators with the radial dependence r~ 3 and 5(r)/r. It was shown in Ref.|| 
that logarithmically divergent contributions are mutually cancelled. This cancellation 
means that for the states with vanishing angular momentum, the correction we discuss 
splits naturally into two contributions - those of large and small distances - each gaining 
its value in its own scale. To calculate the former, one can use again the nonrelativistic 
approach, whereas the short distance contribution, residing in the Compton wavelength 
order scale, calls for a true relativistic approach. 

The closed expression for the first recoil correction to an energy of the relativistic 
electron moving in the Coulomb field is outlined in Sec. 2. This expression is used in 
Sec. 3 for the evaluation of the long-distance contribution. It proves that the relativistic 
approach is more efficient even at the atomic scale. The contribution of short distances is 
found in Sec. 4 employing the Feynman gauge. Sec. 5 is devoted to checks of the obtained 
results. The long-distance contribution is recalculated there using the nonrelativistic 
approach, while the short- distance one is found in the Coulomb gauge. Finally, in Sec. 6, 
we give the numerical values for the energy shifts and compare results of the present 
work with those obtained earlier in [fj, |J. 

Throughout the paper the relativistic units h = c = 1 are used. Since we do not 
discuss radiative corrections, Z is also set equal to unity. 
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2 Methods of Calculation 



One of the perturbation schemes we use at the present work starts from the Schrddinger 
equation in the Coulomb field, when both particles are considered as nonrelativistic in the 
zeroth approximation. To account for relativistic effects, a kind of the operator product 
expansion is built by calculation of scattering amplitudes for free relativistic particles. 
Thus arising effective operators are then the subject for the ordinary perturbation theory 
of the nonrelativistic quantum mechanics. This approach is rather well suited for long- 
distance contributions, which are due to effective operators saturated by nonrelativistic 
region and having non-local kernels. Formerly it was used in calculations of i) logarithmic 
in a corrections to the spectrum of the two-body system ||, ii) the order ma 6 corrections 
to the positronium P levels ||, and iii) the order m 2 a 6 /M corrections to the hydrogen 
P levels [Q. Unfortunately, this approach becomes very tedious being applied to short- 
distance contributions, when effective operators with local kernels are represented by a 
number of diagrams. 

An alternative approach deals with relativistic light particle (electron) moving in the 
field generated by the slow heavy one (nucleus). In the zeroth approximation, the heavy 
particle holds still being a source of the Coulomb field. Wavefunction of the system 
reduces to that of the light particle satisfying the Dirac equation. To first order in the 
heavy particle's inverse mass, the perturbation operator coincides with its nonrelativistic 
Hamiltonian: 



2M 1 ; 

Here P is the operator of a nucleus' momentum. The vector potential A acts at the 
nucleus' site. 

Unfortunately, one cannot calculate an energy shift induced by the perturbation ([!]) 
straightforwardly, i.e. taking merely its average. In fact, the operator ([!]) depends on 
the nucleus' dynamical variables while the argument of the unperturbed wavefunction 
is a position or momentum of the electron. To overcome this difficulty we use the gauge 
invariance of observables in the quantum electrodynamics 0. Being reexpressed in terms 
of electron's variables, the average value of (|XJ) should be retained gauge invariant. The 
new form of the average is now nearly evident: 

&E rec = -l-(^ Up- 3) G e+u (p- 3)) . (2) 



MJ 2m 

Here p is the electron momentum operator, D is the integral operator describing the 
transverse quantum exchange. It has the kerneQ 

ATiactk _ _ k(ak) 
k 2 — u 2 ' k k 2 



what follows we will write often a kernel rather than an appropriate operator for the sake of 
brevity. 
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In (fj), G is the Green's function for the Dirac equation in the Coulomb field, the average 
is taken over a Dirac-Coulomb eigenstate with an energy E. Actually, the "seagull" part 
of (|2|), that of the second order in D, emerges naturally as a counterpart of the A 2 
term from ([]]) when we take the expectation value of this term over fluctuations of the 
electromagnetic field. All the other terms provide the invariance of (|2|) with respect to 
the gauge transformation, 

ip— > exp(icj)(r))ip , D— >D + i\p, <f)]. (3) 

One can easily convince oneself that just the same result can be obtained from the 
formula (11) of @ with the help of the Dirac equation. 

The first attempt to obtain the relativistic expression for the recoil correction to 
hydrogen energy levels was undertaken in Ref.||. Complete expressions for various 
contributions in the Coulomb gauge were originally derived in the framework of the 
quasipotential approach in Ref. [|]. The sum of those contributions can be convinced to 
reduce to the right-hand side of (Q). 



3 Long-Distance Contribution 



Present section is devoted to the calculation of the contribution to the energy shift 
which is saturated by the atomic scale and can thus be called the long-distance one. To 
check the results, two approaches described above have been applied in parallel. Here 
we describe in detail how the second, relativistic, approach works. The procedure of 
comparison with the results of the more cumbersome nonrelativistic approach will be 
postponed until the Section 5. 



3.1 Pure Coulomb Contribution 



In the relativistic approach the pure Coulomb contribution, 

A£c = -h I Is <^+^ = jm < ?(a + - ■ (4) 

can naturally be represented as the sum of two terms ||, 

AE C = (J—)- — (pA_p). (5) 
6 \2M / M y 1 

Here A + and A_ are the projection operators to sets of positive- and negative-energy 
Dirac-Coulomb eigenstates respectively. With the aid of the Dirac equation the mean 
value of p 2 /2M can readily be reexpressed in the following form |§ : 

I p 2 \ m 2 — E 2 m 2 I ( E \ a a 2 \ 

+ 77i7 2 - + -F • ( 6 ) 



2M 2M 2M \ \m r r 
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As for the second term in (|5|), responsible for virtual transitions into negative-energy 
states, the simple analysis shows that it doesn't contribute to the order of interest at the 
atomic scale. Actually, the trivial power counting on the right-hand side of the obvious 
inequality, 



(pA_p) | < 



1 ([p,c]A4p,c]) 



(7) 



4m 2 

where C is the Coulomb potential, shows that at the atomic scale, the product of commu- 
tators is already of the sixth order in a, so that the projector and the wavefunctions can 
sufficiently be replaced by their nonrelativistic counterparts. Since there is no negative- 
energy states in the nonrelativistic approximation, the atomic scale contribution to the 
initial average also vanishes in the order we consider. 



3.2 Magnetic Contribution 

After performing the integration over u, the expression for the single transverse, or 
magnetic, contribution, 

turns into 

Ap a \m){m\ ^ \m)(m\ \ 4na k \ 

AEM = _ M Y I? k + E m — E ~ ? E-E m + k ) IT ) ' (9) 

where denotes the sum over discrete levels plus the integral over positive-energy 
part of the continuous spectrum, while J2- stands for the integral over negative-energy 
continuum. 

For the transverse photon momenta in the atomic region, k ~ ma, one can expand 
the first term in @ to the power series in the ratio (E — E m )/k. To zeroth order (in the 
approximation of the instant exchange), we have: 




The sum of the first term and (§) has very simple form ||, 

m 2 — E 2 m 2 a 2 

2M ~ 2MN 2 ' ( ' 

where the standard notations for the Dirac-Coulomb problem are used, 

iV = y(t + n r ) 2 + a 2 , 7 = \J k? — a 2 , 

n r is the radial quantum number, k = — 1 — al. Notice that ( |TT1) reduces to the lowest 
order result for the states with n T = only. 
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As far as we are seeking only for corrections of the even order in a, the next term of 
the expansion to be considered is 

USr* = -^Re (p^(E m - Ef\m){m\^j = ~^([H, [H,p\] , 

(12) 

where H = dp + f3m + C is the Dirac Hamiltonian in the Coulomb field. This term 
describes the effect of retardation. Implying the corresponding operator by its kernel, 
we have 

[H,p\=a—, [H,[H,fi\ = a ^ 

To the lowest nontrivial order, matrix elements of dPs over positive-energy states can be 
replaced by the appropriate Pauli currents: 



J ret 



a 2 I Airk' 2p'k' - k' 2 Att 2p k + iff x k ' 
M\~k^ 2^ "fc 1 2^ 



(13) 



Here k' = q — k, q = p' — p, while p and p' are the arguments of the wavefunction and 
its conjugated respectively. Being converted to the spatial representation, the average 
above equals 

« 2 / -1 7? + 2ffl\ . , 

Strictly speaking, in (0), the integral over k has an infrared divergent part. It is 
omitted from fll4l) since the photon momenta k ~ ma 2 contribute to the previous order 
correction. To make sure that this is correct, one can regularize the divergency supplying 
the photon with a mass A such that ma 2 <C A <C ma. The term proportional to 1/A and 
omitted from fll4|) can be easily checked to cancel the respective term in the difference 
between @ and the expression obtained from (Qj by the replacement k^\Jk 2 + A 2 . On 
the other hand, just this difference determines the low-energy contribution to the order 
ma 5 /M correction. 

The last of contributions due to the single transverse exchange is generated by virtual 
transitions into negative-energy states and is covered by the last terms in @ and (|10D. 
The inequality similar to (0) shows that the nonrelativistic expansion of the last term in 
(§) starts with the seventh power of a. As for the negative-energy contribution to (|i~0D, 
to the lowest nontrivial order it reduces to 

a 2 Uira„ .jftfTrgV a 2 I Aira^k 4ndk \ 
2mM\ k' 2 ' k 2 / Am 2 M\ k' 2 k 2 / ' 1 ' 

yielding in the spatial representation 

AEm - = -^m\-±-^-)- (16) 

Taken over S states, this average is logarithmically divergent at small distances (linear 
divergencies cancel each other). An ultraviolet divergency will be discussed later for 
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the total long-distance contribution to the S level shift. In fact, due to the gauge 
dependence of an individual contribution (e. g., (|T6|) ), its divergent part alone has no 
physical meaning. 



3.3 Seagull Contribution 



Again, taking the integral over u in the expression for the double transverse, or seagull, 
contribution, 

US* = ~^f^ (DG E+ uD) , (17) 



MJ 2ni 
we obtain 

* E a 2 I 47ra fc / \m){m\ / E m - E \ 4na k \ 

S ~2M\ k' + (E m — E + k')(E m — E + k) V k' + k J k " '/ ' 

(18) 

where the ellipsis stands for the negative-energy part which differs from the positive- 
energy one by the overall sign and signs before k and k! . One can easily check that the 
linear in k/2m terms in the expansion of the negative-energy part cancel each other. But 
just these terms at the atomic scale could produce the energy correction of the necessary 
order. Hence it remains to consider the positive-energy part explicitly written in fllBp. 
In the leading nonrelativistic approximation, 

a 2 Una*, 4ira k \ 

we again replace matrix elements of a's over positive-energy states by the Pauli currents, 

a 2 I Att 2p' k , +iaxk' 4ir 2p k + iff x k \ 
AEs+ ~ 2M \V 2 2^ P 2^ / ' (20) 

and perform the Fourier transformation to obtain 

a 2 / ^1 ^ 1 3p + 2al\ , . 

A£s+ = 4^M 2P ^ P+ ^ P-)- (21) 



3.4 Total Long-Distance Contribution 



Summing up fll4|) , ([Tq) and (pT|) we arrive at 



a 2 I P 4vr(5(f)\ , . 

' 2 ZI + — —)■ (22) 



4m 2 M \ r 4 
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Expanding (|TT| ) to the power series in a 2 and evaluating the average in (|22| ) we obtain 
the long-distance contribution for a state with nonzero I: 



m 2 a 6 f 1 6k 
A£/> = -r^<!:^ + ^^ — — — + 




k|(4k 2 - 1)(2« + 3) 8nft 2 
k 2 (k+1) \ 1 ] 



1)(2k + 3)/ 2n 3 



(23) 



For I — 1 it reproduces the result of the paper 0. As previously mentioned, the 
effective operators of the order under consideration contain singularities insufficient to 
compensate the vanishing of a wavefunction with nonzero I at r— >0. That is why ( |23| ) 
is the total sought-for correction to Z > levels. 



For S states the first term in fl2"2]) evidently vanishes due to the angular momen- 
tum operator / annihilating their wavefunctions. It is interesting to note that a naive 
generalization of the result for states with nonzero I to S ones leads to the error - the 
vanishing of the angular average is compensated by the linear divergency of the radial 
one. 

As for the second term, which formally contains the linear divergency, it is just a 
remnant of the short-distance contribution to the previous (fifth in a) order correction. 
To make sure that nothing is lost in the sixth order, let us regularize the ultraviolet 
divergency by subtraction of the potential generated by the massive transverse exchange 
(the photon mass A ^> ma), from the potential of the ordinary transverse exchange 
entering (|T5]): 

k' 2 ~* k' 2 k' 2 + X 2 ' [ ' 

By going to the spatial representation we obtain the regularized version of the singular 
operator: 

— -> — An5(r). 

r 3 

Being averaged, it gives the energy correction of the order mXa 5 /M. The latter should 
be cancelled by the linear in A term in the expansion of the short-distance contribution to 
the order m 2 a 5 /M correction, calculated with the massive propagator of the transverse 
quantum (actually the expansion parameter is X/m <C 1). Along with linear in X/m 
correction, one could expect the correction linear in a = ma / m. However the expansion 
parameter at large distances is (p/m) 2 ~ a 2 so that the operator we discuss does not 
contribute to the order of interest. On the other hand, a linear in a correction to a local 
(oc S(r)) operator can arise as an ordinary radiative one. In this case the correction is 
completely saturated by small distances. The next section is devoted to the calculation 
of such corrections. 

So, in S states the long-distance contribution is exhausted by the m 2 a 6 /M term from 
the expansion of (JTTJ) , 

ld m 2 a 6 fl 3 2 1 \ 
AEl=0 ~ 2Mrf U + 4^ " ^ + • (25) 
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It vanishes in the ground state only. 



4 Short-Distance Contribution 



As we have mentioned in the Introduction, the long-distance contribution for S states 
is supplied by a short-distance one residing at scales of the Compton wavelength order. 
Since two contributions are well separated, each of them is gauge invariant, so that 
evaluating the short-distance one we can use a different, more appropriate gauge. Rather 
naturally the mostly convenient gauge is the Feynman one. The main formula @ can be 
rewritten in this gauge by application of the Dirac equation or directly from the eq.(ll) 
Qf§: 

a 2 f duj I Ait / k'\ „ / A-\ 4tt \ . . 

AEtot = -mI^i \k" + \ 2 -^ [ a + u) Ge+uj [ a -u) F + p-^/ • (26) 

Momenta of the photons are assumed to flow both from the nucleus to the electron. The 
photon mass A is introduced to establish control over infrared divergences reminiscent 
of lower-order and long-distance contributions. Those divergencies arise in the process 
of the approximate evaluation of the integrals in (|26|) . 

Taking the wavefunctions at the origin and replacing the Green's function by the 
first term of its expansion in the Coulomb field, we have 

m + uj + (3m + ap' 4-7T 



AE 



G 




m + uj + (3m + ap 




*-tAz* 3 ■ (27) 



p 2 -n 2 

Here ip 2 = |-^(0)| 2 , the angle brackets denote integrations over p and p'; q = p' — p; and 



Voj 2 - A 2 , n = V2 



2 



mu + uj 

Contrary to the case of large distances, in the deep relativistic region the opposite order 
of integration is suitable - first over p and p', and then over uj. As for the former, 
it becomes rather trivial after conversion to the spatial representation. Preparatory to 
such the conversion, it is convenient to express all the scalar products containing different 
momenta in terms of their squares. Then some of the denominators can be cancelled. 
At this point we can drop those terms which do not contain Q in their denominators. 
In fact, the only scale leaving in such terms is A so that they cannot produce a short- 
distance contribution. In the spatial representation the initial two-loop integral with 
zero external momenta turns into a simple one-dimensional integral over r. The contour 
of the resulting (^-integration encloses the cut between the points —2m and —A in the 
complex plane. After this last integration we obtain 

A£ ° = 1&-^/- ^^= + ^2-S|, (28) 
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where e = X/2m. 



If the Green's function is taken to zeroth order in the interaction, we have to use the 
Dirac equation in order to account for the momentum dependence in the wavefunction: 

. ^ a 3 ib 2 [ duj I An , 
AE^ = 2 — — / ( — (2m + ap' , 



M J 2m \p' 4 v q 
m + uj + f3m + dp 

p 2 -n 2 




(29) 



Using the same procedure, we take — r/2 as the Fourier transform of An/p 4 . The linearly 
divergent constant we thus leave aside is actually proportional to 1/at and contributes 
to the previous order correction. The result of the integration is 

AE, = ^( 1 + i - * r -*—). (30) 



Mm I 2e 2 e 3iTy^e Ji ,/ x (x 2 — 1 



Regulator-dependent terms in (|2"8D and (|5DD arise from the integrals saturated by the 
region of momenta p ~ A and frequency uj ~ A (or ymA) and are thus the remnants 
of the previous orders corrections or of the long-distance contribution. Truly relativistic 
contribution comes from the region p ~ uj ~ m and does not depend on the infrared 
cutoff: 

f4h,2-5W (31) 



Mn 3 V 2, 

It is pertinent to note here that this result is truly short-distance, i. e. it does not contain 
hidden long-distance contributions, which naively could arise due to cancellation of the 
same nonzero powers of A from numerator and denominator - all positive powers of the 
photon mass were dropped out in the process of calculation. On the other hand, an 
emergence of such contributions would be self-contradictory. Actually, if an integral is 
saturated by distances of 1/A order, then at p ~ A, the integrand denominator has at 
least one power of momentum more than the product of the numerator and the measure 
of integration. In other words, any "long-distance" contribution (determined by the 
scale of A) has to contain a positive power of the photon mass in its denominator. 



5 Checking of the Results 
5.1 Long-Distance Contribution 

To be certain that the long-distance contributions are found correctly, all of them were 
rederived in the framework of the nonrelativistic approach which exploits the Schrodinger 
equation as a starting point. For the states with nonzero angular momenta we used the 
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following procedure. All the contributions prove to have the same analytic structure in 
k, namely 



AE 



ma, 

Mn 3 ' 



where 

E 




« -1/2 « + l/2 (k+1/2) 2 k + 1 k + 3/2 / 
6-1/2 \ 



ra/s; 2 I °° k + 1/2 



+ n 2 \ K \ \ k - 1/2 k + 1/2 k + 3/2/ n 3 ' 1 J 

Constants a, b, c, d evaluated in the nonrelativistic approach for individual contributions 
as their asymptotic values at k^oo or residues at corresponding poles were then com- 
pared with the respective results obtained in the relativistic approach. In the process of 
comparison, a number of 'nonrelativistic' contributions breaks down into groups accord- 
ing to the meaning of respective 'relativistic' ones. For example, the retardation part 
of the magnetic contribution flT4j ) comprises three terms in the nonrelativistic approach: 
AEjmc: £^E>mcc an d ^^ret (notations are from ref.0). 

As we mentioned earlier, S states should be treated separately in order to avoid 
fictious contributions arising due to the compensation between vanishing angular aver- 
ages and linearly divergent radial ones. All the ultraviolet divergencies in S states are 
checked to cancel each other. To this end we regularize the effective potentials which are 
too singular at r— s>0 and ensure that the total long-distance contribution for S states is 
independent of the regularization parameter. 



5.2 Short-Distance Contribution 

In order to compare the results of the present work with those of []|], the short-distance 
contribution was calculated using the Coulomb gauge also. A mass of the magnetic 
quantum was used as the infrared regulator. The scheme of calculation is completely 
analogous to those used previously in the case of the Feynman gauge. The short-distance 
contributions are: 

ira 3 tfj 2 1 

Cg = °* = ~M^2' (33) 



= " M^' n W (35) 
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S G = ^(41n2-2), (36) 

S, = ^2 In A. (37 ) 

* Mm 2m v ; 

Here C, M, S denote Coulomb, magnetic and seagull contributions respectively. It is 
easy to check that the sum of these contributions coincides with (EHf). 



6 Conclusion 



Numerically, the correction to the energy equals 2.77 kHz for the ground state and 0.51 
kHz for the 2S state. Being somewhat less than the naive estimate (m 2 a 6 /M « 10.2 kHz) 
it nevertheless is quite comparable with the accuracy of the near future measurements. 
It is also interesting to note that the corrections to 2S and 2P levels (with the radiative- 
recoil correction to 2P level taken into account) are rather close to each other, so 
that the correction to their difference, 0.04 kHz, can be considered as negligibly small 



at the present level of the experimental accuracy fll0 |. 



Let us now set up a correspondence between the results of the present work and 
those of the other papers. The result for I > levels appears to be firmly established 



nj. For S levels our result is contradictory to the recent results of the analytic 



and numerical [111 calculations. 



To elucidate the origin of the disagreement, we consider the correction to the ground 
state energy. It is easy to verify that our short-distance results (|33"1) -(|37D are in one- 
to-one correspondence with respective "high-energy" contributions from [3J]. A similar 
statement is true for long-distance contributions (low- and intermediate-energy ones in 
notations of M), with one exception. The coefficient -2 from Eq.(68) of Ref.|| for the 
intermediate-energy contribution to the retarded exchange by the magnetic quantum, 
differs from our result, -1 (in the same units m 2 a 6 /M), which arises after trivial averag- 
ing in ([14]) over the ground state. Unfortunately, we have not managed to reproduce the 
coefficient -2 starting from Eq.(67) of Ref.(14). Furthermore, several arguments can be 



brought forward, that the result (68), 0] for the retardation contribution looks at least 
suspicious. In particular, the logarithmic divergency in the order m 2 a 6 /M is known to 
appear due to the relativistic corrections to the instant transverse exchange. The result 
of the present work concerning the origin of this logarithmic divergency and the value 
of the corresponding coefficient is contained in ( |T5D and is in complete agreement with 
those of Refs.|| and (T^j. As for the effect of retardation, it gives rise to the finite con- 
tribution only (in accord with (0)). But it follows from the result of Ref. || that just 
the retardation is the source of not only logarithmic, but even the linear divergency at 
small distances, while the long-distance relativistic correction to the instant transverse 
exchange does not contribute at all. 
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